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JZ* ■ Abstract 

ON . 

\ We study the radiative (j) decay into 7r°7r°7 and 7r°r/7 taking into account mech- 

anisms in which there are two sequential vector- vector-pseudoscalar or axial- vector- 
vector-pseudoscalar steps followed by the coupling of a vector meson to the photon, 
considering the final state interaction of the two mesons. There are other mechanisms 
in which two kaons are produced through the same sequential mechanisms or from 
(j> decay into two kaons and then undergo final state interaction leading to the final 
pair of pions or ir°r], this latter mechanism being the leading one. The results of the 
| parameter free theory, together with the theoretical uncertainties, are compared with 

■ the latest experimental results of KLOE at Frascati. 

> : 

^ ■ 1 Introduction 

The radiative decays of the into 7r°7r°7 and 7r°777 have been the subject of intense study 
[H [21 EH HI 03 US HI- O ne of the main reasons for this is the hope that one can get much 
information about the nature of the /o(980) and a (980) resonances from the invariant 
mass distribution of the two mesons. The nature of the scalar meson resonances has 
generated a large debate pj, to which new light has been brought by the claim that these 
resonances are dynamically generated from multiple scattering with the ordinary chiral 
Lagrangians (H1E3E]- 

These two reactions involving the decay of the <j) are special. Indeed, the does not 
decay into two pions because of isospin symmetry. A way to circumvent this handicap is 
to allow the <fi to decay into two charged kaons (with a photon attached to one of them) 
and the two kaons to scatter giving rise to the two pions (or tt°i]). The loop which appears 
diagrammatically is proved to be finite using arguments of gauge invariance [P^ I ITT ^ IH] . The 
radiative decay through this mechanism was studied in [Sj and the results of lowest order 
chiral perturbation theory (xPT) were used to account for the K + K~ — ► 7r°7r° transition. 
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Since the chiral perturbation theory K + K~ — > 7r°7r° amplitude does not account for the 
/o(980), the excitation of this resonance has to be taken in addition, something that has 
been done more recently using a linear sigma- model in [Tlj . 

The work of [Zj leads to the excitation of the /o(980) in the 7r°7r° production, or the 
a (980) in n°7] production in a natural way, since the use of unitarized chiral perturbation 
theory (UxPT), as in Pj, generates automatically those resonances in the meson meson 
scattering amplitudes and one does not have to introduce them by hand. 

The experimental situation has also experienced an impressive progress recently. To 
the already statistically significant experiments at Novosibirsk [1311111111 one has added 
the new, statistically richer, experiments at Frascati [rH|IT?l] which allow one to test models 
beyond the qualitative level. In this sense although the predictions of the work of [7j, using 
UxPT with no free parameters, provided a good agreement with the experimental data 
of |15)I16)I17| . thus settling the dominant mechanism as that coming from chiral kaon loops 
from the <fi — > K + K~ decay, the new and more precise data leave room for finer details 
which we evaluate in this paper. 

One of the issues concerning radiative decays is the relevance of the sequential vector 
meson mechanisms. The anomalous vector- vector-pseudoscalar couplings [2~U) I21|. followed 
by the coupling of the photon to the vector mesons through vector meson dominance, allows 
mechanisms for reactions of the type that we study, through a sequential V — > VP — > PP7 
process. This mechanism is known to provide the u — > 7r°7r°7 radiative width with accuracy 
[2~T] and has been further extended to study p — > 7r°7r°7 and other radiative decays in 
j22H 23j. In particular, it was found that in the latter reaction the sequential vector meson 
mechanism and the loops (in this case charged pion loops) were equally important, and 
the consideration of the two mechanisms produced a width compatible with the recent 
experimental determination of the SND collaboration [21]. Since the — > pn is OZI 
forbidden, the contribution from the sequential vector meson mechanism should be small 
and it was not considered till recently [25)15] . It proceeds via <j) — u mixing, which is indeed 
small, but still noticeable in the radiative <ft decay width. 

In the present paper we also take into account this sequential mechanism, but in ad- 
dition we consider the final state interaction of the two mesons within the framework of 
UxPT. Furthermore, in the line of considering the final state interaction (FSI) of the 
mesons we also allow the sequential vector meson production of kaons, which interact 
among themselves and lead to the final two mesons. This is easily implemented in the 
coupled channel formalism of (§] • The mechanisms involving final state interaction all con- 
tribute to the production of the / (980) or a (980) resonances, since the resonance poles 
appear in the coupled channel formalism in all the elastic or transition matrix elements. 

Another novelty of the present work is the consideration of sequential mechanisms 
involving the exchange of an intermediate axial- vector meson (J PC = 1 ++ or l + ~), both 
producing directly the final meson pair or through the intermediate production of kaons 
which undergo collisions and produce these mesons. 

All the mechanisms considered here contribute moderately, but appreciably, to the <ft 
radiative width. The inclusion of all these mechanisms leads to results compatible with 
the experimental data of Frascati, particularly if theoretical uncertainties are considered, 
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which is something also done in the present work. 

The good agreement with experiment is reached in spite of having in our approach a 
width for the /o(980) very small, of the order of 30 MeV, seemingly in contradiction with 
the 'visual' /o(980) width in the experiment, which looks much larger. The reason for this 
has been recently discussed in and stems from the fact that, due to gauge invariance, 

the amplitude for the process contains as a factor the momentum of the photon, which 
grows fast as we move down to smaller invariant masses of the two pseudoscalars from the 
mass of the /o(980) where the photon momentum is very small. This distorts the shape 
of the resonance, making it appear wider. We shall see that our approach, which respects 
gauge invariance, introduces automatically this factor in the amplitudes. 

We shall see that there is some discrepancy of the theoretical results with the data at 
small two meson invariant masses. We shall discuss this feature, realizing that the results 
resemble very much the raw data, before the analysis is done to subtract the contribution of 
ujit and to correct for the experimental acceptance. Furthermore, some of the assumptions 
made in the analysis of [TH] might be questionable. 

2 The (f> —> 7r°7r°7 decay 

2.1 Kaon loops from — > K + K decay 

The mechanism for radiative decay using the tensor formulation for the vector mesons has 
been discussed in (2H1IZ| and we briefly summarize it here. The diagrams considered are 
depicted in Fig. Q where the loops contain K + K~ . The vertices needed for the diagrams 
are obtained from the chiral Lagrangian 



(i) 



where the notation is defined in |2Sj- We also assume ideal mixing between the <p and the 



'2 1 



7T 1 - vT 8 = 



(2) 



and the matrix V is then given by 

( 72 p % + 75^ 
V 



K*+\ 

uv 



K* 



'*() 



(3) 



L n ° 4-l,,i K 

K*° 6 i 

The amplitude for the radiative decay can be written as T' iu e fl ('y)e u ((f)) = t, and Lorentz 
covariance demands to be of the type 
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Figure 1: Loop diagrams included in the chiral loop contributions. The intermediate states 
in the loops are K + K~ . 



T» v = ag flu + b P^P V + c P»q v + d P v q» + e q^q" 



(4) 



where P, q are the <fi and photon momenta respectively. Gauge invariance {q v T^ u = 0) 
forces b = and d = —a/(P ■ q). If one works in the Coulomb gauge only the g^ v term 
of Eq. (0J) contributes which is calculated from the d term to which only the diagrams 
b), c) of Fig. [T] contribute. Since two powers of momenta correspond to P v q^ the loop 
integrals with the remaining momenta are actually convergent [T^ ITS*]. Furthermore, if 
the K + K~ — > 7r°7T amplitude is separated into the on shell part (setting p 2 = m 2 for the 
internal lines) and the off shell remainder, it was proved in j3U] that this latter part does 
not contribute to the d coefficient. This information is of much value, since it allows to 
factorize the on shell meson meson amplitude outside the loop integral. The amplitude for 
the process is given by 



t = -V2-^e(0) • 6( T ) 



M^GvG(Mj) + q[^-G v ) G(Mj) 



where / = 92.4 MeV and G is the convergent loop function of 



tK+K-,TT°TT° (5) 

given by 



G{M^Mj 



1 f 


- b)I( 


a,b) 


Ml 




M 2 


a = —f-, 


b = 


m 2 



(6) 



where m is the mass of the meson in the loop, Mj the invariant mass of the two mesons 
and I (a, b) accounts for the three meson loop and is given analytically in Eq. (12) of 30 . 
On the other hand, G(Mj) is the ordinary loop function of two meson propagators which 
appears in the study of the meson meson interaction in [H| and which is regularized there 
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Figure 2: Diagrams for the tree level VMD mechanism. 

with a cut off of the order of 1 GeV. In Eq. (J3J) the tK+K-,n n = ^^k^^ * s the transition 
amplitude with the iterated loops implicit in the coupled channels Bethe Salpeter equation 
(BS) obtained in [Qij. The parameters Fy, Gy, for the vector mesons are obtained from 
their decay into e + e~, p + p~ or two mesons. If one works at the tree level one finds 
F^ p) = 153 ± 2 MeV from p -> e+e" or = 143 ± 5 MeV from p -> p+p,~ decay. On 
the other hand, from 7r + 7r - decay at tree level one finds G$ = 69 MeV or if one meson 
loop is considered G^y = 55 MeV [2H]- Similarly for the one finds Fy* ] = 161 ± 2 MeV 
from — > e + e~ or i*^ = 151 ± 6 MeV from — > p + p~ . From — > K + K~ decay at tree 
level we find Gy = 55 MeV but this value would be 5% higher from — > K°K°. It is 
interesting to recall that when the iteration of loops implicit in the BS equation is used to 
evaluate pion and kaon electromagnetic form factors in [31] one finds a universal value for 
Gy = 55 MeV and similarly Fy = 154 MeV. Since here we use a similar formalism to the 
one used in [3T] the use of this latter values of the constants according to the approach 
used gives us an indication of the level of uncertainty in these parameters. We thus take 
for the calculations Fy = 156 ± 5 MeV and Gy = 55 ± 5 MeV. 

2.2 Sequential vector meson exchange mechanisms 



Following the lines of [22| l23j in the study of p and to radiative decays and the more recent 
of |25tl6] in the decay, we also include these mechanisms here. They are depicted in 
Fig. 2, where we explicitly assume that the — > p°7r° proceeds via the — uj mixing. 
In order to evaluate these diagrams we use the same Lagrangians as in |5|l2*2"] 

Cyy P = ^e^(d,V u d a VpP) 

C Vl = -4f 2 egA,(QV») (7) 

where <> means SU(3) trace, g = -4.41, G = = 0.016 MeV' 1 , Q = diag{2/3, -1/3, -1/3}, 
V (P) the vector (pseudoscalar) SU (3) matrices and e is taken positive. 

In addition we must use the Lagrangians producing the (p-u mixing. We use the for- 
malism of [321 

£^ = 9^0^ (8) 
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which means that the diagrams of Fig. El can be evaluated assuming the decay of the u 
(with mass M^) multiplying the amplitude by e (the measure of the tfi-uo mixing) given by 



Ml - M - 



(9) 



Values of 0^ of the order of 20000 — 29000 MeV 2 are quoted in which are com- 
patible with e = 0.059 ± 0.004 used in |34| 1 which is the value used here. 

The amplitude for the <p(q*) — > ^lipij^ip^lio) decay corresponding to the diagrams 
of Fig. 2 is given by 



-C 



_2y/2 egfG 



2 



e- 



M 2 



P 2 {a} + {b(P)} P' 2 {a} + {b(P')} 



+ 



Ml - P 2 -iM p T p (P 2 ) M 2 — P' 2 — iM p T p (P' 2 ) 



where P = P2 + q, P' = pi + q and 



(10) 



{a} = e* ■ e q* ■ q - e* ■ q e ■ q* (11) 
{b(P)} = -e* -e q* ■ P q- P-e- P e* ■ P q* -q + e* -qe- P q* ■ P + e-q* e* ■ P q- P 

with e* and e the polarization vectors of the <p and the photon respectively. 

At this point it is worth mentioning that the theoretical expression for the V — > P7 

decay widths T v ^ Pl = \aC 2 ( m 9 m v ) with Cj the SU(3) coefficients given in Table 
I of [36J obtained from the Lagrangians of Eq. fij). gives slightly different results to the 
experimental values from |3Tj. For this reason we can follow a similar procedure to that 
used for the r\ — > 7r°77 decay in jH^l where the Cj coefficients were normalized so that the 
theoretical V — > P7 decay widths agree with experiment. In the <p — ► reaction this 

procedure results in including in Eq. (fTIHl a normalizing factor C = 0.869 ± 0.014, obtained 
considering the V — ► P7 reactions shown in Table I of 36J. 

2.3 Pion final state interaction in the sequential vector meson 
mechanism 

Since the irn interaction is strong in the region of invariant masses relevant in the present 
reaction we next consider the final state interaction of the pions in the sequential vector 
meson mechanism. 

We must take into account the loop function of Fig. Ek, but on the same footing we 
must also consider those of Fig. Eb andEt, where charged pions are produced and allowed 
to interact to produce the 7r°7r° final state. The thick dot in Fig. El means that one is 
considering the full tttt — ► tttt t-matrix, involving the loop resummation of the BS equation 
of ref. P and not just the lowest order tttt —> rnv amplitude. 

In order to evaluate those diagrams we must calculate the loop function with a p and 
two pion propagators. First let us note that due to isospin symmetry the ujp°7T° coupling 

1 Note that in a different sign for e is claimed. This is actually a misprint and the results of that 
paper are calculated with e > |35| . 



6 



a) b) c) 

Figure 3: VMD diagrams with final state interaction of pions 

is the same as the up + 7i~ or up^7i + . Next, given the structure of the terms in Eqs. 
and (fTT]) we must evaluate the loop integrals 

dAp p»p» 1 \ \ (12) 



(2tt) 4 P 2 - M 2 + ie (q* - P) 2 - m\ + ie (q - P) 2 - m\ + ie 

For simplicity we evaluate these integrals in the reference frame where the two meson 
system has zero momentum. In this frame the and photon trimomentum are the same, 
q. We need the following integrals, which for dimensional reasons we write as 



<1> P P P°D 1 D 2 D 3 ^I (13) 



(2^) 4 

r d 4 P (f 
7 (2tt) 4 d |gl 1 

r d 4 P n 1 ^ 

1 J W7 plpWlD2D3 = Ma + ~¥ h 

where D\, D2, D3 are the three meson propagators of Eq. (|T2*|) . The integrals are evaluated 
in the Appendix. For that purpose the P° integral is evaluated analytically and the d 3 P 
integral is evaluated numerically by means of the same cut off which has been used to 
regularize the two meson loop in the meson meson interaction, i.e., a cut off of the order 
of 1 GeV. 

The amplitude is then written as 

t = -C^-^e* ■ eq {W° -q)+ I b q*° + I q - W° + q)} 2t£ m (14) 
where 

Ml - M 2 n M 2 + M? 

«=^mr • " = M - + "=^mr- (15) 

Note that, since we evaluate the amplitudes in the Coulomb gauge and e is transverse 
to q, e is not modified by a boost connecting the two meson system and the <fi rest frame. 
The same thing happens to e* 



*+ CO, ( 



+ 



,--'k° \k° 

K *° to, <Sf 



,--'k° K° 



K 



*0 (0, c 



Figure 4: VMD diagrams with final state interaction of kaons 

Thus the invariant t matrix of Eq. (j!4j) is evaluated with the polarizations in the rest 
frame, like the amplitude in sections 12.11 and 12.21 The function multiplying e* • e, which is 
thus invariant, is evaluated in the rest frame of the two mesons, as we have done. 

In Eq. (I14|) we have also taken into account that 

< ttV + 7T+7T- + Tr-TT+ltlvrV >= 2^=°^ (16) 

where is the / = tctt scattering amplitude in the unitary normalization of the states 
used in Note that we have used the prescription in which the inr scattering amplitude 
is factorized on shell in the loops and this is justified in |Hj for the loops implicit in the BS 
equation, in [30J for the loops of Fig. [H as we mentioned before, and in [3B] for looping 
diagrams of the type of those we are now discussing. Note also that we take only the s- 
wave part of the — > Tin scattering amplitude, the p-wave part for 7r°7r° being forbidden. 
Anticipating results, we should mention that in the case of 7r°r) in the final state we shall 
also use only the s-wave scattering amplitude since the p-wave part is zero at lowest order 
xPT and negligible if higher orders are considered (32] . 

2.4 Kaon loops in the sequential vector meson mechanisms 

Next we consider the diagrams analogous to those in Fig. 3 but with kaons in the interme- 
diate states. 

Note that unlike in sections 12.21 and 12.31 the <pVP vertices are now not OZI forbidden. 
They come from the Lagrangian of Eq. ((7|). All the four <pKK* vertices in Fig. 4 have 
the same strength and taking into account the wy, 07 couplings coming from Eq. (|7J) we 
obtain for the set of diagrams of Fig. 0] 

t = -C^I^e* ■ eq {l a (2q*° - q) + I b q*° + I q - h(q*° + q)} ■ (17) 



S 



where the Ij integrals are now evaluated using K, K and K* in the three meson loop. In 
Eq. [T7| we have used that 

< K + K~ +K-K+ + K°K° + K^K \t\n 7i° >= -^=t T =4 



V3 



KK,nn 



with t^^thel = KK — > nir unitarized transition amplitude evaluated with the states 
with the unitary normalization of |5] . Unlike in the case of pion loops where only the u is 
attached to the photon, now we can have the u and the <j) without violating the OZI rule. 
The two terms in the last factor of Eq. (|17j) account for the uj and (f) meson respectively. 



2.5 Sequential axial vector meson mechanisms 

Since the mass of the <fi is around 250 MeV higher than the p mass, and we are considering 
sequential vector meson mechanisms with p or K* exchange, we should pay attention to 
the analogous mechanisms involving vector mesons with a similar mass difference with the 
cf) on the upper side and these are the axial and vector mesons with J PC = l + ~ or 1 ++ 
(see Table 1). Therefore, the bi or ai axial vector mesons and the Kib, Kia strange axial 
vector mesons will play the role of the p or the K* in former diagrams. 





/ = 1 


1 = 


I = 1/2 


1+- 


&i(1235) 


/ii(1170), /ii(1380) 




1++ 


ai(1260) 


/i(1285), /i(1420) 


K 1A 



Table 1: Octets of axial- vector mesons. 

Because of the C parity of the states, the Lagrangians for the axial-vector-vector- 
pseudoscalar couplings have the structure of < B{V, P} > for the b\ octet and < A[V, P] > 
for the octet of the a\ [HI], where the <> means SU(3) trace. In the last expressions V 
and P are the usual vector and pseudoscalar SU (3) matrices respectively and 



B 



A = 



( tM 11 ™) + T2 b ° 
^/i(1285) + 



V2 



■1A 







Kts \ 


/ii(1170) - 




K° 1B 


K° 

IB 




/ii(1380) j 






Kf A \ 


/i(1285) - 




n lA 


n lA 




/i(1420)/ 



(19) 



In addition one has to consider an approximate 50 % mixture of the Kib and K 1A states 
to give the physical ^(1270) and ATi(1400) states [1211111110]. 
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We have modified the original Lagrangian of 140.. to treat the vector fields in the tensor 
formalism of 29.. This formalism has the advantage that without basically changing the 
rates of the A — > VP decays, one deduces the coupling of the a\ to ttj using vector meson 
dominance through a\ — > up — ► tc-/, with an amplitude which is gauge invariant and which 
is in agreement with the chiral structure of [2H] for the a\ — > P7 couplings and with the 
experiment. Details are given elsewhere in [4"3*] . 

We hence use the Lagrangians [Hi] 

£ BVP = D<B^{V^,P}> (20) 
C AVP = tF <A^[V^,P] > 

where the i factor in front of the F is needed in order Cavp to be hermitian. 
In Eq. (J2UJ) the fields W^ u = A pu , B pu are normalized such that 

<0\W flu \W-P,e>=-^- [P„ e v (W) - P v e^W)] (21) 

In addition the propagators with the tensor fields are defined as [2E| 



<0\T{W^W pa }\0>=iV^ pa = (22) 
g pp g va (M w - P 2 ) + g pp P v P a - g pa P v P p - (fx <-> v) 



M w 2 



M w -P 2 -ie 

The physical i^i(1270) and Ki(1400), with a mixture around 45 degrees 2 found in 
E1HH1IIHI , can be expressed, in terms of the 1 = 1/2 members of the 1 + ~(1 ++ ) octets, 
K 1B (K 1A ), as 

K 1 (1270) = ^=(K 1B -tK 1A ) 

tfi(1400) = 4=(#ib + iK 1A ) (23) 

With the values for D = -1000 ± 120 MeV and F = 1550 ± 150 MeV, very similar to 
those found in j4"2*ll41U40| . we are able to describe all the A — > VP decays plus the radiative 
decays of the ai — > ttj |4"Hj . 

Once again the <fi sequential decay at tree level through b\ exchange is OZI violating and 
proceeds via — u mixing. We have evaluated this contribution and found it negligible, 
thus we do not further discuss it. Similarly the loops involving pions are equally negligible. 
On the other hand ai exchange is not allowed by C and G parity. 

2 It is worth mentioning that in [421143] two more possible solutions for the mixing angle around 30 and 
60 degrees were found. This uncertainty will be taken into account in the evaluation of the error band in 
our final results. 
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Figure 5: Diagrams for the sequential mechanisms involving the K\ axial- vector mesons. 

2.6 Kaon loops from sequential axial vector meson mechanisms 

The relevant mechanisms involving axial-vectors are those in Fig. |S] in which K, K are 
created and through scattering lead to the final 7r 7r° state. These are not OZI forbidden 
and have a nonnegligible contribution. 

Since we are using the tensor formulation for the vector mesons this forces us to use the 
tensor coupling of the photon to the vector mesons obtained from the Fy term of Eq. |T|): 



Cv, = -e^\vV° u (d»A»-d»A») 
with Ay = 1, |, — *y for V = p, u, <p respectively. 



(24) 



Now the amplitude with ^(1270) as intermediate state is given, following the lines of 
sections 12.31 and 12.41 as 



AeF v 



q* I + I a {2q-q^)+I b q-h{q + q 



V3 



,i=g 

KK, 



,*() 



[D 2 - F 2 
QV2M 2 



V2{D + Ff 



*0' 



6M| 

+ (q- q*°)G KR (Mj) 



(25) 



K!=Ki(1270) 



The amplitude for the diagram with intermediate Xi(1400) is the same as Eq. (|25|) but 
replacing F by — F and using the 1^(1400) propagator in the evaluation of the Ii integrals. 

The term proportional to G K ^{Mi) at the end of the former expression comes from the 
M\ — P 2 of Eq. (|2~2"Jl which cancels the intermediate vector meson propagator, thus leaving 
the loop with just two pseudoscalar propagators, which is the same one that appears in 
the meson meson scattering. 
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Figure 6: Diagrams for the tree level VMD mechanism for the — > ir ^ decay. 

3 The (f> —> 7r°?77 decay 

After the discussion of the former points the consideration of the <p 7i"°^7 decay requires 
only minimal technical details which we write below. 

3.1 Kaon loops from — > KK decay. 

The diagrams to be considered are exactly the same as those in Fig. changing the final 
7T°7r by 7i°r]. Technically all that we do is to substitute in Eq. (J3J) the amplitude t K + K - j7T o n o 
by the tx+K-^r] or equivalently — -j^t 1 ^ n0jl in the nomenclature of 

3.2 Sequential vector meson contribution 

We have now the diagrams 

The amplitude of this contribution is very similar to Eq. (|10jl . We find 

■ _ r ~ 4 egftf f P» + {b(P)} P*{a} + {b(P')} ] 

^3 M. 2 [M 2 -? 2 - iM p T p (P*) Ml - P' 2 - iM^T^P'Z) J 1 ] 



3.3 ttt] final state interaction of the sequential vector meson mech- 
anism 

The results for the amplitude corresponding to FSI in diagram is easily obtained from 
Eq. ()14|). We multiply the expression by y/E, change the masses and widths of the corre- 
sponding intermediate particles in the loops and change 2t£^ r7r by £™iL. The amplitude 
corresponding to the FSI in diagram is done in the same way but we multiply by — Vd 
instead of a/6 to account for the negative sign in Eq. (f2T)J) . 

3.4 Kaon loops from sequential vector meson mechanisms 

The diagrams involved are the same as in Fig. H] but instead of u and <fi coupled to the 
photon we have now only the p and the final 7r°7r° is substituted by ir°T]. Altogether one 
must change to in the u term of Eq. (fTTj) and substitute n7T by — 3\f\t I KR -k V - 
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3.5 Axial vector meson exchange mechanisms 

The tree level with 61, a\ exchange and their octet partners are also OZI forbidden and 
negligible like in the case of 7r°7r° final state. However, the term with kaon loops equivalent 
to Fig. El but with 7r°rj in the final state and a p coupled to the photon instead of u, <ft are 
OZI allowed and we take them into account. 

Once again the amplitude is easily obtained from that of the u part of Eq. by 
changing to and t 1 ^ wn to — ^^t^k-irq and the same considerations regarding the 
contribution of the i^i(1270) and Ki(1400) intermediate states. 



4 Results 



4.1 Differential cross section 

Using the transition amplitudes described in the previous sections, we can calculate the 
differential decay widths of the <fi meson as, 



dwr^L du M^^ 0(1 - cos 



(27) 



where Mj is the invariant mass of the final two mesons, m! is the pion mass for the 7T7T7 
decay and the rj mass for the irrff decay and q is the photon momentum in the initial <ft 
rest frame. O is the step function and 7r o 7 accounts for the angle between the n° and the 
photon and it is given by 



cos 0_o 



7T U 7 



1 

2pq 



{M4 - u{p) - q) 2 - rri 



rl 



2 2 

P ~Q 



where p and u(p) are the 7r° momentum and energy in the initial rest frame. A symmetry 
factor 1/2 must be implemented in Eq. (|27|) in the case of 7r°7r° in the final state. 

The spin sum and average of the transition amplitudes, J2 \ t\ 2 , can be expressed using 
the contravariant tensor F l i as: 



EE 1*1 



where the tensor expression i™ of the transition amplitude t is defined as 

t = F iJ e .(y)e j (j). 



(29) 



(30) 



4.2 Results for the 



7r°7r°7 decay 



First of all, we show in Fig. [7| the contribution of the chiral loops (solid line) together with 
the sequential VMD mechanisms at tree level of Fig. El (dotted line) and the VMD with 
final state interaction, Fig. (dashed line). In the dashed-dotted line we have plotted the 
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Figure 7: Several contributions to the 7r°7r° invariant mass distribution for the — > 7r 7r°7 
decay. Solid line: chiral loops of Fig. Dotted line: sequential VMD at tree level of Fig. |21 
Dashed line: sequential VMD including the FSI of Fig. El Dashed-dotted line: Coherent 
sum of all these mechanisms. (Experimental data from |18j). 



contribution of all these mechanisms together. We can see, that the effect of the VMD 
tree level contribution and its unitarization is small compared to the chiral loops from 
cf) — > K + K~ decay. The total contribution is still rather different than the experimental 
values of [TB] . although the inclusion of the sequential vector meson mechanism provides 
some strength for the distributions at low invariant masses where the kaon loops from 
cf) — ► K + K~ are negligible. 

In Fig.lHJwe can see the contribution of the sequential mechanisms involving kaon loops 
(Fig. EJ), separating the (dashed line) and to (dashed-dotted line) contributions, and the 
coherent sum of both amplitudes (solid line). We can see that the sum of both mechanisms 
almost cancels, thus making the contribution of the kaon loops of the sequential VMD 
mechanisms very small, although the individual contribution of the loops with an uj or a 
meson attached to the photon are sizeable when they interfere separately with the chiral 
loops from — > K + K~ decay, as it can be seen in Fig. El Had this accidental cancellation 
not happened the contribution of the kaon loops of the sequential mechanisms would have 
been sizeable. 

In Fig. El we show the contribution of the axial- vector meson resonances. We can see 
the contribution of the mechanisms with a i£i(1270) as intermediate state (dashed line) 
and the one with ^(1400) (solid line) (see Fig. |HJ). The size of these mechanisms by 
themselves is very small but when they interfere with the chiral loops from — > K + K~ 
decay they give an important contribution, as shown in Fig. ^2 Thus, these mechanisms 
with the axial-vector mesons in the intermediate states cannot be neglected. 

In Fig. ^1 we show the contribution of the different relevant mechanisms as they have 
been added to the model. In the dashed line we can see the contribution of the chiral 
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Mj (MeV) 



Figure 8: Contribution of the sequential mechanisms involving kaon loops. Dashed line: 
(j) contribution. Dashed-dotted line: uj contribution. Solid line: Coherent sum of <fi and u) 
contributions. 




Figure 9: Addition of the sequential loops contribution to the chiral loops of Fig. ^ Solid 
line: chiral loops of Fig. ^ Dashed line: chiral loops of Fig. ^ + sequential loops involving 
kaons with (j). Dashed-dotted line: chiral loops of Fig. ^ + sequential loops involving kaons 
with uj. Dotted line: Coherent sum of all these contributions. 



15 



0.7 r 

0.6- 



2 05 " 

I-h 0.4- / \ 

1 0.3- / \ - 

o*0.2- / \ - 

A - / I - 

©0.1- / 1 - 

b i 

•a - ^ x i " 

300 400 500 600 700 800 900 1000 



M, (MeV) 



Figure 10: Contribution of the sequential mechanisms involving axial-vector mesons. Solid 
line: i^i(1270) contribution. Dashed line: i^i(1400) contribution. 




Figure 11: Addition of the sequential loops involving axial- vector mesons to the chiral 
loops of Fig. [T] Solid line: chiral loops of Fig. [TJ Dashed line: chiral loops of Fig. [T] + 
axial-vector meson contribution. 
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Figure 12: Different contributions to the two pion invariant mass distributions of the 
<f) — > 7r°7r°7 decay: Dashed line: chiral loops of Fig. ^ Dashed-dotted line: chiral loops 
of Fig. Q + sequential VMD and its final state interaction. Solid line: idem plus the 
contribution of the mechanisms involving axial- vector mesons, (full model). 



c 

0.869 ±0.014 


e 

0.059 ± 0.004 


G v (MeV) 
55 ±5 


F v (MeV) 
156 ±5 


U (MeV) 
92.4 ± 3% 


A (MeV) 
1000 ±50 


D (MeV) 
-1000 ± 120 


F (MeV) 
1550 ± 150 



Table 2: Parameters which uncertainties are relevant in the error analysis. The f n and A 
are the /„■ constant and cutoff of the momentum integral respectively in the loops involved 
in the unitarized meson-meson rescattering. 

loops. In the dashed-dotted line we add the contribution of the sequential VMD and 
its final state interaction. In the solid line we add the contribution of the loops of the 
sequential mechanisms with axial-vector mesons in the intermediate states. This latter 
line represents the full model. 

Up to now, all the curves shown in the figures have been calculated using the central 
values of the parameters without considering the uncertainties in their values. In Fig.lT3*lwc 
show the final result but including an evaluation of the error band due to the uncertainties 
in the parameters of the model. This error band has been calculated implementing a 
Monte Carlo gaussian sampling of the parameters within their experimental errors. The 
parameters of the model which uncertainties are relevant in the error analysis are shown 
in Table H 

The errors in f n and A assumed in the calculations have been chosen such that the 
quality of the fit to the im phase shifts along the lines of is still acceptable within 
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Figure 13: Final results for the 7r°7r° invariant mass distribution for the — > 7r 7r°7 decay 
with the theoretical error band. Experimental data from |18| . 



experimental errors. 

The parameter with the larger contribution to the error band turns out to be the 
Gy since the largest contribution, chiral kaon loops form <p — > K + K~ decay, is roughly 
proportional to Gy (up to the term with q(^- — Gy) in Eq. (j3j) which would be zero within 
some vector meson dominance hypotheses [2H] and is small with our set of parameters). 

The total width and branching ratio obtained in the present work are 



r^c V o 7 = 520 ± 150 eV (31) 
BR((j) -> 7r°7r°7) = (1.2 ± 0.3) x 10~ 4 



to be compared with the experimental values 



BR exp ((f) -> 7r°7r°7) = (1.22 ± 0.10 ± 0.06) x 10" 4 |TB] 
BR exp ((p -► 7r°7r°7) = (0.92 ± 0.08 ± 0.06) x 10" 4 [To] 
BR exp ((f) 7r°7r°7) = (1.09 ± 0.03 ± 0.05) x 10" 4 [IB] 

In Fig. we can see that our results, considering the error band, fairly agree with the 
experimental data except in the region around 500 MeV. The reason of this discrepancy 
will be further discussed in Section |S] 

4.3 Results for the <j) — > tt t]j decay 

In Fig. El we show the contribution of the chiral loops from <fi K + K~ decay (dashed- 
dotted line) together with the sequential VMD at tree level (dotted line) and the loops of 
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Figure 14: Several contributions to the 7T r] invariant mass distribution of the — > 7r°r]7 
decay. Dashed-dotted line: chiral loops of Fig. [T] Dotted line: sequential VMD at tree 
level. Dashed line: loops of the sequential VMD involving kaons. Solid line: coherent sum 
of all these mechanisms. (Experimental data from 

the sequential VMD involving kaons (dashed line). In the solid line we show the coherent 
sum of all these contributions. 

In this case the VMD mechanism involving ni] creation followed by ni] — > tii] FSI 
is negligible and, therefore, is not shown in the figure. We can see that the kaon loops 
of the sequential VMD mechanism give a very important contribution, in contrast to the 
cf) — ► 7r°7r°7 case where the important contribution of the VMD sequential loops was the one 
involving pions. Recall that the reason for the small contribution of kaon loops in the VMD 
mechanism for — > 7r°7r°7 was the accidental cancellation of the and to contributions (for 
the and to attached to the photon). However, in this case there is only the p attached 
to the photon and thus the contribution of this term is large, as was also the case in 
— > 7r°7r°7 for the individual contributions. 

In Fig. El we show the contribution of the mechanisms with intermediate 1^(1270) 
(solid line) and 1^(1400) (dashed line). 

These contributions by themselves are very small but when they interfere with the 
chiral loops give a sizeable contribution, (from dashed to solid line in Fig. IT5j) . 

In Fig.^lwe show the different relevant contributions when they are added one by one. 

The dotted line represents the chiral loops by themselves. In the dashed line we have 
added the contribution of the sequential VMD mechanisms together with the loops in- 
volving the FSI of the mesons. Finally, in solid line we have add the contribution of 
the axial-vector contributions. This latter curve represents the full model. In this figure 
we can see the importance of the loop mechanism involving kaons and the nonnegligible 
contribution of the mechanisms with axial-vector mesons in intermediate states. 

Again, in Fig. we have plotted the full model performing the theoretical error analy- 
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Figure 15: Contribution of the sequential mechanisms involving axial-vector mesons. Solid 
line: Xi(1270) contribution. Dashed line: i^i(1400) contribution. 




Mj (MeV) 

Figure 16: Addition of the sequential loops involving axial-vector mesons to the chiral 
loops of Fig. [T] Solid line: chiral loops of Fig. [TJ Dashed line: chiral loops of Fig. [T] + 
axial-vector meson contribution. 
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Figure 17: Different contributions to the 7r 7] invariant mass distributions of the — > 
7r°^7 decay: Dotted line: chiral loops of Fig. ^ Dashed line: chiral loops of Fig. ^ + 
sequential VMD and its final state interaction. Solid line: idem plus the contribution of 
the mechanisms involving axial- vector mesons, (full model). 

sis 3 . We can see that when these uncertainties are considered we obtain a theoretical band 
in acceptable agreement with the experimental data. 
The total width and branching ratio obtained are 

Y^fi m = 250 ± 80 eV (32) 
BR((p -> ttV) = (0.59 ± 0.19) x 10~ 4 

to be compared with the experimental values 



BR exp (<p -»■ 7r%7) = (0-88 ± 0.14 ± 0.09) x 10" 4 [T7j 
BR exp ((f) -> vrV) = (0-90 ± 0.24 ± 0.10) x 10" 4 HHj 
BR exp ((p -> ttV) = (0-85 ± 0.05 ± 0.06) x 10" 4 HH 

5 Further discussion of the (f> —> 7r°7r°7 results. 

We would like to comment on the strength that we obtain around 500 MeV in the 
invariant mass distribution in the <fi — ► 7r°7r°7 decay which appears in contradiction with 

3 We have also checked that the use of a mixing angle for the strange members of the axial nonets of 
around 30 or 60 degrees |42U43| turns out in decreasing the lower limit of the error band in around 5% and 
10% for the <j> — > 7r°7r°7 and <j> — > n°r]^ decays respectively. 
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Mj (MeV) 

Figure 18: Final results for the n°ri invariant mass distribution for the <fi — > tt -!]^ decay 
with the theoretical error band. 

the experimental analysis. As we saw it comes from accumulation of the novel mechanisms 
which we have discussed in our paper. Such mechanisms are not considered in other 
theoretical papers which find a good agreement with the data. For instance in jB] a good 
reproduction of the data is obtained using a linear sigma model with <r(500) and /o(980) 
included in the calculation and some free parameter to adjust to the data. In [27J an 
attempt to make a model independent calculation is done by parametrizing amplitudes 
compatible with analyticity and gauge invariance and, through a fit to data, a set of 
parameters leading to agreement with the data is obtained. In the same approach 
for the chiral kaon loops from (f> — > K + K~ decay as in the present paper is followed, but 
in addition a contact term in (p'yK !^ is added and, fitting the strength of this term to 
the data, a good agreement with the experimental results is also found. Our philosophy 
has been different and we have not fitted any parameter to the <fi radiative decay data 
but simple have considered the different mechanisms that can sizeable contribute to the 
process. An acceptable agreement with the data is found in the region of the / (980), 
which is the most important issue concerning this reaction. This is not trivial a priori in 
view of the very small width of the dynamically generated /o(980) (around 30 MeV) that 
one obtains in the model of that we use here and the large 'visual' width of the / (980) 
peak in the present experiment. Part of the reason for the agreement comes from the 
q factor (photon momentum) in the amplitude, as requirement by gauge invariance and 
emphasize in |2l?l l2T]. which gives more weight to the amplitude as we move down in the 
TT7C invariant mass from the upper limit (where q — 0). However, as seen in our results, the 
inclusion of the new mechanisms and their interference with the dominant one, particularly 
the contribution of the axial-vector meson exchange mechanisms, also contributes to the 
widening of the distribution around the / (980) peak. 

Although the agreement with data at low masses is not very good, we must point out 



22 



two sources of uncertainty in the experimental spectrum. First, the results in the low and 
intermediate mass region largely depend on the background subtraction dominated by the 
non-resonant ujit process. The size of this process is difficult to obtain because it has 
a strong background itself, mostly from the — > /o7 process, as it is discussed in |3E] . 
There, its magnitude has been obtained in a model dependent way assuming some a priori 
spectrum for the — > f j process jl^l- In fact, before the subtraction, the raw data 
resemble much more our calculated spectrum, (see fig. 4 from [IB])) and we could think of 
a slightly smaller ujtt background. 

Additionally, there is some uncertainty in the way the data are corrected to account 
for the experimental efficiency. This is done in ^BJ by dividing the observed spectrum by 
the effect of applying the experimental efficiency on some theoretical distribution. This 
unfolding procedure depends on the theoretical model used, which we think at low 7r°7r° 
masses is at least incomplete. In fact, with the unfolding method used, the zero value of 
the spectrum obtained with the theoretical model of ^H] implies unavoidably a zero value 
for the corrected experimental results. A reanalysis to the light of the present discussion 
would be most welcome. 

6 Conclusions 

We have studied the radiative decay of the into two neutral pseudoscalars, tc h , Tr r], 
and a photon. We have taken an approach to study the process addressing mechanisms 
proved to be relevant in the study of the radiative decays of the p and the u. These include 
the kaon loops from the —>■ K + K~ decay and mechanisms of sequential vector meson 
steps with vertices of the type VVP, together with vector meson dominance to couple a 
vector meson to the photon. 

In addition to these mechanisms studied before we have added the final state interaction 
in the sequential vector meson mechanisms, allowing for rescattering of the 7r°7r° or 7i°r] 
states. Simultaneously we also allowed kaons to be produced in the sequential vector meson 
processes, followed by the interaction of the kaons to lead to the final 7r°7r° or TT°r] states. 
These latter mechanisms have the advantage over the tree level sequential vector meson 
mechanism that they are not OZI forbidden. We found these kaon loop contributions to 
be individually large, although they become small in the — > 7r°7r°7 decay due to an 
accidental cancellation of two terms, and they are very important in the tt°^7 decay. 

Another addition in the present work is the contribution of sequential axial-vector 
mechanisms involving two sequential AVP steps followed by the coupling of a vector meson 
to the photon, and even more important the production of KK through these mechanisms 
followed by final state interaction of the kaons to give 7r°7r° or ir°rj in the final state. 

The final state interaction of pairs of mesons involving diagrams with loops has been 
done using techniques of unitarized chiral perturbation theory which allow one to study 
meson meson interaction in coupled channels up to 1.2 GeV and hence are particularly 
suitable for the present reaction. 

The main contribution to the radiative decay is the loop mechanisms involving the 
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interaction of K + K~ coming from decay, which by itself reproduces qualitatively the 
experiment as it was already shown in [JJ . However, the new mechanisms studied here are 
by no means negligible. Altogether we find a good agreement with experiment, particularly 
when the study of the theoretical uncertainties is done. These theoretical uncertainties, 
stemming from errors in the magnitudes used as input to determine the parameters of the 
theory, are of the same order of magnitude as the experimental ones. Apart from these 
uncertainties we should stress that there is no freedom in the theory used. There are no free 
parameters which have been adjusted to the — > 7r°7r°7 and — > n ^ data studied here. 
The agreement with the experimental data both for 7r°7r° and ir°r] channel is remarkable 
and less than trivial when one realizes that a change of sign in one of the mechanisms leads 
to opposite interference effects and to results in not so good agreement with the data. The 
agreement is more remarkable, or surprising, when one realizes that this agreement with 
the data is reached for the 7r°7r°7 case, in spite of having an / (980) resonance that has a 
width much smaller than the apparent width that one could guess from the — > 7r°7r°7 
experiment, and which one would get fitting the data with models less elaborate than the 
present one. Some explanation for this can be already seen in [2*^ 12*7] which remarked that 
the requirement of gauge invariance in the amplitude introduces the factor of the photon 
momentum which tends to widen the distribution of the invariant mass around the /o(980) 
peak. Our approach is manifestly gauge invariant and one can show explicitly in all the 
terms that this factor appears. Our approach comes to stress the point of view of |2T^l2*7] 
that the apparent width of the /o(980) in the experimental data of this reaction can not 
be used as an direct measure of the actual / (980) width. 

The main conclusion of the work is that a theoretical approach to the problem of the 
— > 7r°7r°7 and — > Tr ^ decays is possible by using the chiral unitary approach used 
successfully to reproduce the im data up to 1.2 GeV and many other physical processes 
As we could see along the work, the / (980) or a (980) have not been explicitly coupled 
to the 0, meaning that there is no direct — > /07, a o1 coupling. This is in accordance 
with the philosophy of the chiral unitary approach about these two resonances which are 
dynamically generated by the meson meson interaction. This has as a consequence that 
in processes where these resonances are produced one should not include them explicitly 
in the approach but they should appear naturally as soon as the meson are allowed to 
interact. In this sense the present work comes to reinforce the claim that the /o(980) and 
ao (980) are dynamically generated resonances and that the meson meson chiral Lagrangian 
at lowest order contains the driving forces that make possible this generation through the 
multiple scattering of the mesons in coupled channels. 

Acknowledgments 

Two of us, J.E.P. and L.R., acknowledge support from the Ministerio de Education, Cultura 
y Deporte. This work is partly supported by DGICYT contract number BFM2000-1326, 
and the E.U. EURIDICE network contract no. HPRN-CT-2002-00311. 



24 



_ P-q* 



M, 



P-q 



Figure 19: Three meson loop. 

Appendix: Integrals involving three meson propaga- 
tors 

For the evaluation of the loop of Fig. HH1 we need to calculate the integrals of Eq. (JT3j) . 
From Eq. (fTHj) it is easy to obtain that 



la 



d 4 P 

d 4 p 

i r d A P 



2 J (2tt) 4 
% r d A P 



P°\P\ cos eDtDzDs 
P 2 sin 2 8DiD 2 D 3 



2 J (2tt) 4 
i 



P 2 (3 cos 2 9 - l)DiD 2 D 3 



i n _ 



and D 3 - {P _ q)2 _ m , +ie . 



(33) 



where D\ 

Defining tu v = y i P 2 + My, cu 2 = \j {P — q) 2 + m| and ui = y 1 {P — g*) + m 2 , the 
analytical structure in the complex P° plane of the integrals of Eqs. (}3*3*j) is depicted in 
Fig. EDI 

Applying the Residues Theorem with the path shown in Fig. |2Dl the evaluation of the 
P° integral gives 



A; 



_ 1 rdcos6d\P\ \P\ 2 

2 J (2vr) 2 uj^ 2 ( g *o + UJl+ Uv ) ( g *o _ UJl _ UJv + ) ( g *o _ q o + UJl + U j 2 _ ie j 



1 



(q° + lu 2 + LU V )(q° -lu 2 -lu v + i^^)(q° - q*° + u) x + uj 2 - ie) 
for i = {0,1, a, 6}, with s v = (q*° - Uxf - P 2 , s' v = (q° - u 2 f - P 2 and 



(34) 
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Figure 20: Poles and path for the evaluation of the P° integral. 
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